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Our Letter [1] shows that there is an equivalence between the Dicke model and the infinitely 
range XX (Lipkin) model, at the level of the effective Hamiltonian describing the qubit reduced 
subsystem. The Comment by Brankov et al. asserts that our use of a reduced Hamiltonian 
in place of the original Dicke Hamiltonian, is "wrong" as compared to the results which they 
previously reported. 

Brankov et al.'s criticism is misguided - in short, they are wrong to use the word "wrong" 
since we make no statement, anywhere in our Letter, to claim that our effective (i.e. reduced) 
qubit Hamiltonian is exact at finite temperatures. Yet it is on this basis, and this basis alone, 
that Brankov et al. use the word "wrong". There is nothing wrong with our Letter. As a 
side note, we emphasize that their own earlier work is effectively irrelevant to the core results 
of our paper. In addition, there are several other confused statements in the Comment by 
Brankov et al. which, for the sake of completeness, we point out below. 

First, the goal of our Letter was to study the entanglement features of the qubit subsystem, 
not the entire qubit-radiation system. We focused on the quantum phase transition which, 
by definition, occurs at zero temperature. In particular, we looked at the finite-size scaling of 
two-qubit entanglement and contrasted it with previously published numerical results for the 
full Dicke Hamiltonian [2]. Our main finding was that even at the level of the lowest-order 
approximation (i.e. summing special parts of all cumulants) our reduced effective Hamiltonian 
at zero temperature yields a quantum phase transition at the same critical coupling constant as 
the full Dicke Hamiltonian. Even more importantly, our reduced effective Hamiltonian shows 
the same critical finite-size scaling exponents for two-qubit entanglement (i.e. concurrence). 
In addition, we found a universal scaling function for the behavior of the concurrence. 

Second, at no point in our Letter do we dispute Brankov et al.'s claim that our temperature- 
dependent Hamiltonian (Eq.(8) in [1]) is not exact. Obviously it is not. Instead, and as stated, 
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it is the result of a lowest-order approximation. In short, we never presented it as an exact 
result for non-zero temperatures. However, what is remarkable about our result is that at 
zero temperature it does coincide perfectly with the exact one, and also with the finite-size 
quantities which can be obtained using it. This is acknowledged by Brankov et al. in their 
Comment, hence they do agree with our results. We emphasize that although the cumulant 
approach which we used [3] can be worked out to any order, the expressions of higher order 
terms turn out to be analytically quite complex - hence our finding that lowest-orders can be 
quantitatively accurate at zero temperature, is very important. 

Finally, our work is the first to have shown explicitly that the critical scaling exponents 
of two-qubit entanglement (i.e. concurrence) are indeed identical for two different physical 
models (i.e. Dicke and Lipkin models). In addition, we showed that the quantum phase 
transitions of both models belong to the same universality class. Brankov et al. do not 
dispute this finding, nor there is any proof in their Comment to invalidate our basic result. 
The question raised concerning the critical behaviour of finite-size systems remains an open 
problem of course - but this was beyond the scope of our Letter. 
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